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Abstract. 

The problem of the relationship between entanglement and two-qubit systems in 
which it is embedded is central to the quantum information theory. This paper suggests 
that the concurrence hierarchy as an entanglement measure provides an alternative 
view of how to think about this problem. We consider mixed states of two qubits 
and obtain an exact solution of the time-dependent master equation that describes 
the evolution of two two-level qubits (or atoms) within a perfect cavity for the case of 
multiphoton transition. We consider the situation for which the field may start from 
a binomial state. Employing this solution, the significant features of the entanglement 
when a second qubit is weakly coupled to the field and becomes entangled with the 
first qubit, is investigated. We also describe the response of the atomic system as it 
varies between the Rabi oscillations and the collapse-revival mode and investigate the 
atomic inversion and the Q-function. We identify and numerically demonstrate the 
region of parameters where significantly large entanglement can be obtained. Most 
interestingly, it is shown that features of the entanglement is influenced significantly 
when the multi-photon process is involved. Finally, we obtain illustrative examples 
of some novel aspects of this system and show how the off-resonant case can sensitize 
entanglement to the role of initial state setting. 
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1. Overview 

As a lot of researchers suggest [1-5], entanglement is one of the most mysterious 
aspects of quantum physics. Although it is by now possible to verify the predictions 
of entanglement theory in a variety of experiments [2], there remains a considerable 
gap between the formal definition of entanglement and the observable effects that are 
associated with this property. A major thrust of current research is to find a quantitative 
measure of entanglement for general states. For the experiments in the newest fields of 
physics, quantum computing, quantum communication, and quantum cryptography [6- 
13] the quantitative analysis of the multi-qubit or ion is of substantial interest [14]. The 
dipole-dipole interaction between two atoms can be understood through the exchange 
of virtual photons and depends on the transition dipole moment of the levels involved. 
It can be characterized by complex coupling constants, or by their real and imaginary 
parts, where the former affect decay constants and the latter lead to level shifts [15,16]. 
There is an inherent interest in analytical and non-perturbative solutions of multi-atom 
interacting with the cavity field problems, all the more considering quantum systems 
with more than one qubit. One examples of such kind is the system of two two- level 
qubits in an electromagnetic field [14-26]. Entanglement of identical qubits is a property 
dependent on which single-qubit basis is chosen, as any operation should act on each 
identical qubit in the same way. Indeed, individual qubits are excitations of a quantum 
field, and the single-qubit basis defines which set of qubits are used in representing the 
many-qubit state [27]. The analysis of entanglement sharing of the two-atom Tavis- 
Cummings model has been discussed in Ref. [28]. Recently, interest has mounted in 
exploring the quantum system composed of two qubits interacting with a thermal field 
[29]. We have addressed a general two-qubit system in a recent paper [23] in which an 
analytical expression for the temporal evolution of the Pancharatnam phase when the 
field starts from vacuum state is given. 

This has motivated us to criticize the conception of entanglement of a two-qubit 
system in the context of the mixed quantum states. The main contribution of this 
paper is to synthesize conceptual insights that already exist to push forward a more 
coherent view of how concurrence as an entanglement measure of a two-qubit system 
might make progress in understanding the two-qubit entanglement. To be more precise, 
we assume that two two-level atoms (two qubits) share a bipartite system, taking into 
account the multi-photon transition. Another principal aim is to elucidate the extent 
to which mixed entangled states can affect the entanglement. The emphasis being 
put on the investigation of the entanglement in a more general situation in which the 
two atoms (qubits) share a mixed state, rather than a pure state and we propose to 
use the concurrence hierarchy as a measurement of entanglement [30-32]. The issue 
of attributing objective properties to the constituents of a quantum system composed 
of two qubits, does not turn out to be a straightforward generalization of the just 
analyzed case involving distinguishable qubits, and the problem of entanglement has 
to be reconsidered. Entangled mixed states may arise when one or both qubits of an 
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initially pure entangled state interact, intentionally or inadvertently with other quantum 
degrees of freedom resulting in a non-unitary evolution of the pure state into a mixed 
state. In general it is known that there are also cases when entangled states are mixed 
with other entangled states and where the sum is separable. 

The rest of the paper is organized as follows. We devote section 2 to give a brief 
overview of the binomial states and initial states setting. In section 3, we present 
notations and definitions of the model and its analytical solution to be used in the 
rest of the paper. Section 4 is devoted to consider the atomic inversion as well as 
the quasiprobability distribution. The concurrence as an entanglement measure is 
presented in section 5, followed by a numerical computation in which we shall examine 
the influence of different involved parameters in subsection 5.1. The paper is closes with 
the conclusions outlined in section 6. 



2. Binomial state 



The implementation independence in quantum information theory is guaranteed by 
the use of Hilbert spaces, states and operations between and on them. It is not said, 
what they physically describe in more concrete terms, whether we are dealing with 
spins, polarizations, energy levels, qubit numbers, or whatever you can imagine. The 
features of nonclassical states visualize specific aspects of nonclassicality and do not 
yield a complete characterization of nonclassicality as a phenomenon on its own. In 
this section we give a brief overview of the binomial states [33]. For measuring the 
quantum state of the radiation field, balanced homodyning has become a well established 
method, it directly measures phase sensitive quadrature distributions. Alternatively, 
unbalanced homodyning yields access to phase-space functions. For some systems, such 
as a cavity-field mode and a mode of the quantized motion of a trapped ion, methods 
have been proposed that allow for a direct detection of the characteristic functions of 
the quadratures. The binomial states are finite linear combinations of number states 
[33] 
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|n>, (1) 



which interpolate between some fundamental sates such as number states and coherent 
states, where to is a non-negative integer, rj is areal probability (0 < r\ < 1) and \n) is 
a number state of the radiation field. Such states have been studied in great detail in 
the literature (see e.g. Ref. [34]). The binomial states have the properties 

{| to) rj — > 1 
|0) (2) 
| a) rj — > 0, to — * oo, rjm = a 2 . 

In the theory of open system or the reduction theory, one often considers two subsystem 
A and F represented by Hilbert space. Let &j, (j = A, F) be state spaces (the set of 
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all density operators) and &a ® &f denotes the state space in the composite system. 
Here, we assume that, before entering the cavity, the field initially in a binomial state 
such as 

zu = \rj, m)(r),m\ E & F . (3) 

In pure-state quantum mechanics the state of the system is usually represented by 
a normalized wavefunction, which is a unit vector in a Hilbert space. Entangled 
mixed states may arise when the qubits of an initially pure entangled state interact, 
intentionally or inadvertently, with other quantum degrees of freedom resulting in a 
non-unitary evolution of the pure state into a mixed state. In general it is known 
that there are also cases when entangled states are mixed with other entangled states 
and where the sum is separable. The usual interpretation of mixed states, is that 
their creation involves irreversibly destroying information [35,36]. This has interesting 
consequences concerning entanglement theory, since there, the irreversibility is often 
associated with the fact that one is dealing with mixed states. We assume that, the 
two two-level qubits initially prepared in the mixed states enter the cavity whose single 
mode under consideration, prior to the interaction with the qubits, is in the binomial 
state zu. Thus, the atomic density matrix is of the diagonal form 

p? (0) = cos 2 ^ \a) ( ei \ +sm 2 9 i \g i ){g i \ E & Al , (4) 

where p^(0) is the atomic density matrix for the %— qubit. The initial state of the two 
qubits system can be written as 

p a (0)=p?(0)®p«(0)G6 A . (5) 

It is well known that mixed states can be realized by an ensemble of pure states in 
an infinite number of ways. The determination of the separability of a state and 
the determination of its entanglement have in common that a particular realization 
of a state has to be found such that some property holds for all pure states in that 
realization. In order to find this optimal realization, it is of considerable interest to 
have a mathematically elegant way of generating all possible realizations of a state [31]. 

3. Model for pair of qubits 

We consider a mechanism through which a system of two qubits can be entangled in 
a cavity field. We assume that the qubits are modeled by two-level systems having 
multiphoton transition, which is a micromaser system [21]. The micromaser is an 
experimental realization of the idealized system of a two-level qubit interacting with 
a second quantized single- mode of the electromagnetic field [22]. In this section we 
consider a theoretical model which differs from the standard micromaser set-up in that 
instead of a single qubit we have assumed a pair of qubits interacting with a single mode 
of the cavity field (see figure 1). In the dipole and rotating wave approximation, one 
can write (h — 1) 
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Figure 1. Schematic of the multiphoton transition in the interaction of a pair of two 
two-level qubits with a single-mode radiation field. In the shown process we denote by 
\ e j)(\9j)) the j— qubit's upper (lower) level and {\i}} are the virtual levels. 

+ + + /^V 2) }e(fc) + Ti^Wa* 

+ 7 1 *^ ) ^ + 7 2 ^ ) a fc + 7 2 *a tfc ^ 2 e ) . (6) 

We denote by the atomic operators for the z— qubit. a and a) are field operators 
corresponding to annihilation and creation of photons in the cavity mode. We denote 
by 7j the coupling constant for the i— qubit, uj is the field frequency and u>i is the atomic 
frequency for the i— qubit. The parameter Q(k) is defined such that 

If k is larger than unity, in a consistent physical treatment one should take into account 
field-induced level shifts which are proportional to the number of photons in the field 
mode. Here we denote by pf* and f$ the intensity-dependent Stark shifts to the two 
levels of the i— qubit, that are due to the virtual transitions to the intermediate level. 

We write the Hamiltonian H into two mutually commuting parts H = Hq + Hi n , 
where [H , H in ] = 0, 

H = uj (ata + s<£ + - s$j - sj£) , 

+ A [fi£> + s£) - s£) - s£>] + 7l ^a* + ^5g)&t* 

+ 72^^ + 72^^. (8) 
The detuning parameter is given by A = uji — ku). The continuous map £^ describing 
the time evolution between the qubits and the field is defined by the unitary evolution 
operator generated by H in such that 

e; : e A ^e A ®& F , 

£; P =U t (p a (0)®w)U:. (9) 
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The interaction Hamiltonian, in this case, leads to an exactly solvable time evolution 
operator. Resuming our analysis, the time evolution operator can be written as 

U t = exp I ~ I H(t')dt' I . (10) 



o 



h 

A solution of equation can be written as 

4 4mm 



i=l z=l n=0 1=0 



where 



n iz (n,v,t) = ^(n)^\n)^(v)^(v)nUn,v,0)e- u ^' x ^, (12) 

l,j,s,p=l 

= |ei,e 2 ,n- A;), |* 2 ) = |#i, e 2 ,n), \^ 3 ) = \e 1 ,g 2 ,n), 
1*4) = \gi,92,n + k). 

We denote by A m the eigenvalues of the Hamiltonian H in and 3^(n) is the i— element 
of the I— eigenvector. The coefficients Qi Z (n,v,0) specify the initial conditions for the 
field and atomic states. The eigenvalues A/ are to be found from a fourth-order scalar 
equation, the roots of which may be easily written in closed form for two identical atoms, 
arbitrary detuning in the absence of the Stark shifts or for two nonidentical qubits at 
exact resonance. For a general multiphoton interaction in the presence of Stark shifts 
and detuning parameter, this problem can be treated numerically. If we consider the 
dispersive approximation in which A >> /i n 7, (where 7 = 72 A is the dipole coupling 
constant) and the second qubit is weakly coupled to the field, an explicit expression for 
the final state £ t *p can be easily obtained [37,38]. In that case, it is straightforward to 
obtain explicit expressions for Q iz (n,v,t), namely, 

Q u = A^Af'bn-kbi-k cos 2 9 1 cos 2 9 2 + B^B^h sin 2 8 1 cos 2 9 2 , 
fiia = AlBfb^k-k cos 2 6 1 cos 2 6 2 - S*Af6 n 6, sin 2 6 1 cos 2 6 2 , 

n 22 = A^A\b n bi sin 2 e, cos 2 e 2 + s*s**6 n _ fc 6,_ fc cos 2 e 1 cos 2 e 2 , 

^33 = A^AH^k cos 2 6 1 sin 2 6 2 + B^Bll^n+kh+k sin 2 X sin 2 2 , 
fi 34 = - A^. k Bf +k b n bi cos 2 9 X sin 2 9 2 - B^AllkK+kk+k sin 2 9 X sin 2 9 2 , 
K.+k-A-1+kbn-kbi-k 



fi 44 = A^A^bn-kk-k sin 2 0i sin 2 9 2 + B* n+k B%. k b n bi cos 2 X sin 2 9 2 , (13) 



where 



1 f 72 / 72 

A n = - exp (-it{g n - fi n }) < 1 + h exp (-2i^ n t) 1 - 



2 1 v r J/ [ 2/i n 1 v ™ ' V 2/i. 

Tl / 77/! 

5 n = - — W J exp (-it{£„ - /i n }) {1 - exp (— 2i/x rt t)} , (14) 



A^n = \/(7274)+ w! ' (?„ = A + 72 U+ 2 
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Qij — Qj i7 r2i3 — fii4 — f2 23 — f2 24 — 0. 

A downside of analyzing more complex atomic system is that analytic expressions 
for the final state functions and, consequently, for the matrix elements are not always 
available. Therefore, the greatest benefit of equation (10) which represents an analytical 
solution of the final state of the system for this general model, is that it directly yields 
a method for actually calculating any property related to the system. 

4. Atomic inversion and field properties 

We mainly devote the present section to consider the atomic inversion from which the 
phenomenon of collapse and revival can be observed [18], and see how it is affected in 
the present model. The population inversion expressions for each qubit can be written 

as 

(<r«(*)> = ~7>{|ei) ( ei \ - \ 9i ) ( 9i \} Pi (t), i = 1,2, (15) 

where pi(t) is the reduced atomic density matrix of the i— qubit which can be obtained 
by tracing out the field variables i.e., 

Pi(t) = Trj {Sip) . (16) 

In this case, the total atomic population inversion is given by 

<*(*)> = \ «<x (1) (*)> + ^ (2) (t)». (17) 




10 20 30 40 50 60 10 20 30 40 50 60 

Scaled time Scaled time 



Figure 2. Time history of the total atomic inversion (<i(t)) as a function of the scaled 
time 7ii. Calculations assume that k = 1,02 = ir/4,9i = 0, m = 70,72/71 = 0.2, the 
detuning parameter A has zero value, and for different values of r\ where (a) rj = 0.2 
and (b) 7] = 0.7. 

Figure 2 is a time history of the total atomic inversion (cr(t)) according to equation 
(I17|) against the scaled time 7 1 t. In this figure we assume that the one-photon transition 
k = 1, the second atom starts from a maximum entangled state (#2 = f) an d the first 
atom starts from a pure state 9\ = 0, while initial field is a binomial state with m = 70, 
and for different values of the parameter rj, where 77 = 0.2 for figure 2a and rj = 0.7 for 
figure 2b. The detuning parameter A has a zero value and — = 0.2. In this case and with 
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a small value of rj (say rj = 0.2) the total atomic inversion exhibits well-known collapses 
and revivals. However, any change of the binomial parameter rj leads to changing in 
the atomic inversion and consequence, increasing the parameter rj, leads to elongating 
the revival time while the atomic inversion oscillates around the same value (see figure 
2b). In general when rj increases, the number of isolated revivals decreases at the same 
period of time. A natural next question would be: Given a quantum state of which 




123456 123456 
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Figure 3. The evolution of the atomic inversion (<x(t)) as a function of the scaled time 
7it. Calculations assume that k = 2 (the two-photon processes), 62 = 7r/4, #1 = 0, 
m = 70, 72/71 = 0.2, the detuning parameter A has zero value, and for different values 
of rj where (a) r = 1, 77 = 0.2 and (b) r = 0.7, rj = 0.2. 
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Figure 4. The evolution of the atomic inversion (cr(t)) as a function of the scaled 
time jit. Calculations assume that k = 1,02 = tt/4, 8\ = m = 70,72/71 = 0.2, the 
detuning parameter A has zero value, and for different values of 77 where (a) r\ = 0.2 
and (b) i] = 0.7. 

one knows that it is entangled, how can (<j(t)) be affected by a multiphoton transition? 
One may envision, for example, the situation that a two-photon process is involved (i.e. 
k = 2). In figure 3, we consider this case in which the Stark shift will be taken into 
account (we set k = 2 and the other parameters are the same as in figure 2). We see 
that the total atomic inversion has rapid oscillations with a periodical collapses and 
revivals. This discussion has clearly demonstrated that the general behavior of the total 
atomic inversion, when the two-photon transition is involved, is remarkably affected by 
changing the number of quanta and (<j(t)) is almost periodic. This periodicity follows a 
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consequence of the fact that the generalized Rabi frequency in the two-photon transition 
is proportional to n rather than to ^Jn which is the case in the single-photon process. In 

the presence of the Stark shifts namely r = \J (3^ / (3^ , we note that the Rabi frequency 
as well as the temporal width of the oscillations packets decrease (see figure 3b). The 
most obvious difference, is that the population inversion collapses to a nonzero value 
with the inclusion of Stark shifts but to a zero value without it. This fact highlights an 
important difference between the two situations. 
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Figure 5. Contour plot of the field reduced density matrix quasiprobability 
distribution function against X = Re.(f3) and Y = Im.((3). Calculations assume 
that k = 1, #2 = 7r/4, #i = 0, m = 70, 72/71 = 0.2, and the detuning parameter A has 
zero value, where (a) 71 1 = r/ = 0.2 and (b) 71 1 — r\ — 0.7. 

In what follows we shall choose the parameters to show the first atom starts also 
from a mixed state, i.e Q\ — |. In this case the evolution of the atomic inversion 
still exhibits Rabi oscillations, showing collapses and revivals. The oscillations are 
independent of 9i, but 9\ strikingly affects the amplitude of the inversion (see figure 4). 
As soon as we take 0\ effect into consideration it is easy to realize that the amplitude 
of the oscillations is decreased. Furthermore if we take larger values of the binomial 
parameter 77, the revival time is elongated and the atomic inversion oscillates around 
zero. 

At this end and after discussing a particular aspect of the atomic inversion in the 
multi-photon two two-level qubits model with the initial field prepared in the binomial 
state, which clearly exhibits the collapses-revivals phenomenon and provide us with 
information about the discrete nature of the quantum qubit-cavity field eigenvalues, 
we are going to focus our attention on the representation of the field in phase space 
which provides some aspects of the field dynamics. Perhaps the most convenient 
quasiprobability to be used in this kind of problem is the Q— function. The first step 
to be taken is the calculation of the reduced density operator for the field PF(t), and 
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then we get the Q— function as 

Q(x,y,t) = -((\p F (t)\C), (18) 

7T 

where \() is a coherent state with amplitude ( = x + iy. The Q— function is not only 
a convenient tool to calculate expectation values of anti-normally ordered products of 
operators, but also gives us a new insight into the mechanism of interaction in the 
model under consideration. Figure 5 shows a contour plot of the reduced density matrix 
Q— function quasiprobability for different values of the binomial parameter r\. When 77 
increases, the Q— function bifurcates into two blobs rotating in the complex coherent 
state parameter plane in the clockwise and counterclockwise directions with the same 
speed. The collision of the blobs corresponds to the revival. This can be clearly seen if 
we compare the numerical calculations for inversion and Q— function. 

5. Concurrence 

The characterization and classification of entanglement in quantum mechanics is one 
of the cornerstones of the emerging field of quantum information theory. Although an 
entangled two-qubit state S^p is not equal to the product Sfpi and £^p2 of the two 
single-qubit states contained in it, it may very well be a convex sum of such products. 
In general it is known that microscopic entangled states are found that to be very stable, 
for example electron-sharing in atomic bonding and two-qubit entangled photon states 
generated by parametric down conversion. Entanglement as one of the most nonclassical 
features of quantum mechanics is usually arisen from quantum correlations between 
separated subsystems which can not be created by local actions on each subsystem. 

As more and more experimental realizations of entanglement sources become 
available, it is necessary to develop different methods of measuring the entanglement 
produced by different sources [39]. An ensemble is specified by a set of pairs 
{(pi, \i>i))}iLi, consisting of N state vectors \ipi) and associated statistical weights pi, 
and, N is called the cardinality of the ensemble. The concurrence of a bipartite state 
(i.e., a state over the bi-partite Hilbert space (5,4 <S> &f) is defined by an almost magic 
formula [30,31] 

Cs tp it) = max {0, Ai - A 2 - A 3 - A 4 } . (19) 

We denote by X { the square roots of the eigenvalues of (£ t *p) x i^tP) m descending order, 
where 

{ZtP) = Ki ® <t w ) (S*p)* (a yi <g> a m ) . (20) 

where a Vi is the Pauli matrix. The importance of this measure follows from the direct 
connection between concurrence and entanglement of formation. It has been shown that 
[30] the entanglement of formation of an arbitrary state E^p is related to the concurrence 
Ce* p (t) by a function 

E F (C etP {t)) =7r + (t)log7T + (t)+7r-(t)log7r-(t), (21) 
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where 

^(t) = \(l + yfl ~ C\ tp (t)) . (22) 

The entanglement of formation is monotonically increasing with respect to the increasing 
concurrence. If we could write the general bipartite pure state as [40,31] 

d— 1 oo 

W*))= $>*>>*) K, J, n>, (23) 

i,j=0 n=0 

in this case, the concurrence can be calculated as 
Ce. p (t) = {2[l-Tr{p\(t))]Y 

2 y 

\pik(n,t)p jm (n,t) - pi m (n,t)pjk(n,t)\ > . (24) 




The concurrence Cs* p (t) as a measure of the degree of entanglement ensures the scale 
between and 1 and monotonously increases as entanglement grows. A note of caution 
about how to interpret the state of a physical system in terms of quantum entanglement 
may be in place here. The previous standard definitions of quantum entanglement 
tacitly assume that, every state in the bipartite or multipartite Hilbert space is in 
principle available as a physical state and local as well as global quantum operations, 
measurements and unitary transformations, can be performed on the Hilbert space. 



5.1. Results 

We now discuss applications of the above equation to specific situations. As stated 
above, an important situation is that, when Ce* p (t) = the two qubits are separable 
and Cs* p if) = 1 indicates maximum entanglement between the two qubits. In our 
numerical examples we have used physical parameters from some recent experiments 
[22, 41], but extrapolated the qubit transit time t to rather large values of "fit. It 
is, of course, an experimental challenge to obtain a one-qubit source and atomic life- 
times of the atomic states involved such that these large values of 71 1 can be reached. 
An interesting question is whether or not the entanglement is affected by the different 
parameters of the present system with the initial state in which one of the qubits is 
prepared in its excited state and the other in the mixed state. In particular, the mixed 
state parameter 9 2 , the scaled time "fit, and the parameters from the initial state of 
the field (to, 77). A numeric evaluation of the concurrence as an entanglement measure 
leads to the plot in figures 6-10. We now pause to touch on certain concurrence features 
when i] rs (i.e., the coherent state), the mean photon number rjm = 20 for figure 
6a and rjm = 10 for figure 6b. The maximum value of the entanglement increased 
as the mean photon number is decreased (see figure 6), but the entanglement vanish 
as the time goes on for large values of the mean photon number this is not the case 
when n takes small values (see figure 6b). It is interesting to note that, the maximum 
entanglement is achieved when 8 2 = r! f, n = ±1,±3,±5,.. while Cs* P {t) ~ 0, for 
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Figure 6. Plot of the concurrence Cg* p (t) as a function of the scaled time 71 1 and 
62- The field initially in the binomial states with 77 = 0, k = 1,A = 0,72/71 = 0.2 , 
and for different values of the mean photon number rjm — n , where (a) n = 20 and 
(b) n = 8. 

$2 = 22 f E , m = 0, ±1, ±2, ... Also, the entanglement shows symmetry around 82 = 0. Put 
differently, with different values of rj, such as r/=0.7, one will have a very small amount 
of entanglement at the initial period of time only and this amount disappear when the 
time goes on (see figure 7a). Indeed, the comparison of plots figure 6a and figure 7a 
where 7] = and rj = 0.7, respectively demonstrates that the entanglement in both cases 
has somewhat similar behavior corresponding to different values of 62- In accord with 
the initial conditions 77 = 0.7 and m = 70, figure 7b shows more oscillations at the same 
period of time while the entanglement survive in this case longer. Also, small amount 
of entanglement is repeated several times with the time development. At the period 
2 < 7it < 3, the entanglement shows small oscillations throughout the manipulations 
as a result of increasing the binomial parameter rj = 0.9. The comparison between 
figure 7a and figure 7b shows the obvious effects of the binomial parameter rj in the 
dynamics of the system. Next, we will analyze the influence of dispersive approximation 
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Figure 7. Plot of the concurrence Cs* p (t) as a function of the scaled time 71 1 and 82- 
The field initially in the binomial states with n — 20, k = 1, A = 0, m = 7O72/71 = 0.2, 
and for different values of r\ where (a) r\ = 0.7 and (b) r\ = 0.9. 

on the appearance or disappearance phenomenon of the entanglement previously found. 
Dispersive effects can be conveniently incorporated by assuming that 72 << 7i, (such as 
— = 0.01). It is remarkable to see that with 77 = 0.7, entanglement is nearly washed out 
for the initial stage of the interaction time. While more oscillations have been observed 
when the time goes on (see figure 8a). The situation is changed for 77 = 0.9 (see 
figure 8b), in this case the maximum entanglement increased further and start earlier. 
The initial period in which the entanglement washed out in figure 8a decreased with 
increasing the parameter 77. Also, the maximum value of the entanglement is increased 
with increasing rj. The number of quanta effect on the entanglement is particularly 
pronounced as this number takes large values. In figure 9a the entanglement is plotted 
as a function of 71 i and 82 for k = 2. However, we note from this plot that the 
maximum amount of the two-qubit entanglement indeed does move closer to the point 
at jit = 0. Again we notice some similarities with other plots in the sense that there is 
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Figure 8. The same as figure 2 but 72/71 = 0.01. 

no entanglement for 2 = ™ and maximum value at 9 2 = ^f, (see figure 9a). From our 
further calculations (which were not presented here) it is clear that we can get the same 
amount of entanglement using the present measure (concurrence) and the negativity 
as a measure of entanglement, in agreement with our previous result [42]. Also, the 
qubits and radiation subsystems exhibit alternating sets of collapses and revivals due 
to the initial mixed states of the qubits and radiation employed here. Let us now 
consider the situation when the qubits are initially both in mixed state and nonzero 
detuning. It is surprising that a nonresonant case (nonzero detuning) can entangle 
two qubits depending on their atomic initial preparation. First of all, we note that the 
entanglement appears only for the initial period of the interaction time. These properties 
show that the role played by the detuning on the entanglement is essential. 

The above results pose two intriguing questions: (i) Why does the entanglement 
become maximum for B2 — ^f, and takes zero value for 82 = 0, (ii) Why does the 
entanglement due to the concurrence behave essentially in a similar way with different 
values of the mixed state parameter whatever values of the other parameters? In what 
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follows, we propose an analytic expressions which can give an answer to these questions. 
We now first analyze the reason why the concurrence does not exceed zero value when 
9 2 = 0. To prove that analytically, one may first consider 9{ = 0, then we have only 
nonzero values of f2 n ,fi 12 and Q 22 and all the other f^- vanish, i.e. the coefficients f^- 
in equation (fT3~|) reduce to 

^11 = A^A^bn^kbi-k, 
^12 = A^B^bn-kbl-k, 

n 22 = 6 n _ fc 6,_fe. (25) 

Using these coefficients with the above definition of the concurrence we easily find that 
Cs* P (t) = 0, this is also in agreement with the numerical calculations. On the other 
hand, if we consider the first qubit in its excited state and the second qubit in a mixed 
state, then equation (fT3|) reduces to, 

n n = A^Af 6 n _ fc 6{_ fc cos 2 9, 
Q 12 = AlAfK^b^k cos 2 9, 
n 22 = B t ri Bj*b n _ k b l ^ k cos 2 e, 
n 33 = A^A^bnk sin 2 9, 
fi 3 4 = - Ai^B^bnh sin 2 9, 

flu = B t n+k B t l * +k b n b l sin 2 9. (26) 
For simplicity we used 9 instead of 9 2 . In this case the concurrence is given by 

I sin 29 v 



Ce tP {t) 



I b 2 n b 2 n _ k A^B^f, — B^A 1 *^ I . (27) 



n=0 

From equation (}2Tj) we see that for 9 = we will get maximum value of Cg* p (t) 
and for 9 = r | L , Cs* p (t) = 0. It is also clearly from this equation, the concurrence 
depend on the amplitude of the initial binomial state. Here, we clarify that it can be 
done by using the different initial state, which is strongly affected by the atomic number 
representation. This naturally leads to the use of occupation numbers of different single- 
qubit basis states in quantifying multi-qubit entanglement even when the number of 
qubits is conserved. The occupation-numbers of different modes have already been used 
in quantum computing [43]. It is important to refer here to the work in Ref. [44] in which 
experimentally a superposition state of the ground state and a non-maximally entangled 
antisymmetric state in two trapped ions has been realized. In the experiment two 
trapped barium ions were sideband cooled to their motional ground states. Transitions 
between the states of the ions were induced by Raman pulses using co-propagating 
lasers. The non-maximally entangled state was used [45] to demonstrate the intrinsic 
difference between quantum and classical information transfers. The difference arises 
from the different ways in which the probabilities occur and is particularly clear in terms 
of entangled states. 

Before we conclude, it is necessary to give a brief discussion on the experimental 
realization of the present model. It was reported that the cavity can have a photon 
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Figure 9. Plot of the concurrence Cg* p (t) as a function of the scaled time 71 1 and 
62- The field initially in the binomial states with n = 20, m — 70, r\ = 0.7, 72/71 = 0.2 
, where (a) k = 2, A = and (b) k = 1, ^ = 10. 

storage time of T = 1 ms (corresponding to Q = 3 x 10 8 ). The radiative time of the 
Rydberg atoms with the principle quantum numbers 49, 50 and 51 is about 2 x 10~ 4 s. 
The coupling constant of the atoms to the cavity field is 2n x 24 kHz [46] . The experiment 

may be described schematically as follows: a stream of pairs of rubidium atoms in 
circular Rydberg states, time and velocity selected, was sent through a resonant cavity 
[47]. Each pair consisted of one atom in a state with principle quantum number 51, 
and the other in a state with principle quantum number 50, the frequency of transition 
between the two states being 51.1 GHz. Also, there was an experiment [47] in which 
the phase of oscillation of each of two entangled qubits is individually controlled and 
can be adjusted to provide the necessary values for testing a Bell inequality. 
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6. Conclusion 

To sum up, in the main part of this paper we have investigated the properties of the 

dynamically emerging entanglement in the multiphoton two two-level qubits. We have 

treated the more general case where initial states of the two qubits can be mixed with a 

binomial state of the field. We have obtained an exact solution of the density operator 

taking into account the dispersive limit, and thus provides insight into the behavior 

of more complicated two-qubit systems. We have investigated the interaction of the 

binomial states (which reduce to number and coherent states in two different limits) with 

atomic systems in the framework of the two two-level qubits, and describe the response 

of the atomic system as it varies between the Rabi oscillations and the collapse-revival 

mode and investigate the total atomic inversion and the quasiprobability distributions. 

The idea of using the concurrence as an entanglement measure offers many attractive 

features. Entanglement is measured via the concurrence, currently used only for an 

arbitrary system of two qubits, but similar analysis can in principle be applied to other 

systems such as a bipartite system with arbitrary dimensions. As we anticipated, this 

system exhibits some novel features in comparison with the single qubit system. We have 

found that some different regimes occur, depending on the actual initial joint product 

state, the number of quanta and the binomial state parameters. It was demonstrated 

that quantum entanglement is stored in the model system considered, and that the 

nature of this entanglement is strongly dependent on the detuning of the atomic levels, 

an effect that may have important consequences in other nonlinear processes. 
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